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We present an exact rewriting of the Mie coefficients describing the scattering of light by a 
spherical core-shell particle which enables their interpretation in terms of an hybridization of the 
two surface modes arising, respectively, at the core-shell and the shell-medium interface. For this 
particular case we thus obtain from the Mie theory-analytically for all multipole orders and hence 
for arbitrarily sized particles-the hybridization scenario, which so far has been employed primarily 
for small particles in the electrostatic approximation. To demonstrate the strength of the rewriting 
approach we also extract the hybridization scenario for a stratified sphere directly from the expansion 
coefficients for the electromagnetic fields. 

PACS numbers: 42.25.Bs, 42.25.Fx 


I. INTRODUCTION 

Ever since the pioneering work by Mie 1 and Deby^, 
the classical optical response of objects plays an im¬ 
portant role in applied science 34 . The applications 
range from the spectroscopy of grains embedded in 
gaseous interstel lar environments 5 to plasmonic devices 
on the nanoscal^ 7 . In particular the latter is a grow¬ 
ing field of research driven by the progress in materi¬ 
als synthesis and processing® which provides plasm onic 
nanostructures of continuously increasing complexity ® 
nanoricc 11 !, nanorings 12 -1, or nanoshelld 1 to name only a 
few. The geometry of these structures provides efficient 
means to tailor their optical response, opening thereby 
new fields of application, for instance, in biomedicine 15 
or photovoltaics-^, which in turn stimulates the design of 
even more intricate structures. 

Different methods have been developed over the 
years to analyze light scattering by composite objects. 
While the approached- 1 7 extending the early analyti¬ 
cal workdI 2 'are usually restricted to spherical particles, 
numerical approaches are now available which do not suf¬ 
fer this limitation 25 2 '. Depending on the discretization, 
the methods are either surface- or volume-based. In par¬ 
ticular the latter are very powerful since they can han¬ 
dle arbitrarily shaped inhomogeneous objects. But like 
the generalizations of Mie’s original approach numerical 
methods provide no physical picture of the interaction of 
light with composite objects. 

An appealing p hysica l picture was first given by Pro- 
dan and coworker^® 29 for metallic nanoshells and later 
by Preston and Signoreli 30 for dielectric core-shell par¬ 
ticles. As pointed out by these authors the optical re¬ 
sponse of composite entities such as nanoshells can be un¬ 
derstood within the surface mode hybridization scenario. 
The essence of the scenario is that the surface modes aris¬ 
ing at the interfaces of these objects interact with each 
other and result in optical resonances not to be found in 
homogeneous objects. Modifying the object’s geometry 
changes the interaction and hence the optical response. 
Prodan and coworkers showed that the interaction can be 
interpreted as a kind of hybridization. In analogy to elec¬ 


tronic states of condensed matter, optical resonances of 
composite objects can thus be labelled bonding and anti¬ 
bonding depending on their symmetry. Both Prodan and 
coworkers 28 and Preston and Signoreli 30 worked out the 
hybridization scenario within the electrostatic approxi¬ 
mation using, respectively, a Lagrangian and an eigen¬ 
value method for the description of the charges induced 
at the interfaces of the object. The two approaches are 
very flexible, can be applied to arbitrary geometry, and 
require not necessarily bulk dielectric functions. They 
can be combined with microscopic models for the optical 
response of the atomic constituents of the objects a nd ar e 
thus well suited for applications in nanoplasmonics 9 10 . 

The interpretation of light scattering by composite ob¬ 
jects in terms of the hybridization scenario is often re¬ 
ferred to as a new conceptual approach? 1 or as a funda¬ 
mentally different way of thinking about plasmonic ef¬ 
fects 32 . As the Mie theory is the exact analytical de¬ 
scription of light scattering by spheres and thus naturally 
includes all observed effects, the question arises whether 
the hybridization scenario can be also found directly in 
the Mie formulae and thus whether the new thinking can 
be united with the old formalism. 

In this paper we show that this is indeed the case. 
For the particular case of a spherical core-shell particle 
the hybridization scenario can be straightforwardly de¬ 
rived fr om th e Mie theory by recasting the Mie coeffi¬ 
cients 3 - 4 17l: ® in a form resembling the diagonal elements 
of a matrix resolvent describing two hybridized energy 
levels. The derivation, valid for all multipole orders and 
thus applicable to arbitrarily sized particles, starts with a 
splitting of the electromagnetic fields inside the shell into 
two partJ^ 3 -31, corresponding to the penetrating fields of 
an homogeneous particle made out of shell material and 
the scattering fields of a cavity filled with core mate¬ 
rial and sitting inside an homogeneous domain of shell 
material. In a second step the expansion coefficients ob¬ 
tained by enforcing the boundary conditions at the core¬ 
shell and shell-medium interface are then rewritten into 
the desired form. The re-organization of the fields inside 
the shell and the re-ordering of the expansion coefficients 
it leads to yields no physical concepts beyond the hy- 
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bridization scenario 28 30 . In fact it was inspired by it. 
But it is comforting to see the scenario emerge directly 
from the formulae of the Mie theory, which are notori¬ 
ously difficult to interpret. We thereby also generalize 
the work of Ruppin 19 ' and Uberoi 20 who showed, again 
only in the electrostatic approximation, that the optical 
response of a core-shell particle arises from two coupled 
subsystems. The hybridization scenario can be also found 
in the Mie formulae for a stratified sphere. Essential is 
again the splitting of the fields into penetrating and scat¬ 
tered parts and the re-organization of the expansion coef¬ 
ficients guided by the structure of the diagonal elements 
of a matrix resolvent describing hybridized energy levels. 

In the next section we rewrite the Mie coefficients 
for a spherical core-shell particle and demonstrate how 
the building blocks of the hybridization scenario, the 
medium-embedded homogeneous shell particle and the 
shell-embedded core cavity, can be extracted from the 


new expressions. For illustration we present in Sect. Ill 


numerical results for a dielectric core-shell particle, the 
type of particle we suggested to employ in a gas discharge 
as an electric probe with optical read-out 84 . In partic¬ 
ular we show, up to the hexapole, data for the splitting 
between the bonding and antibonding resonances as well 
as the electric fields inside the particle. The stratified 
spherical particle is discussed separately in an Appendix. 
Concluding remarks are given in Sec. |IV[ 



Figure 1. Geometry of light scattering by a core-shell particle 
with total radius and core radius n embedded in a medium. 
Refractive indices A r , characterize the core (i = 1), the shell 
(i = 2), and the medium i = m. The incident electromagnetic 
plane wave is described by a Poynting vector «Si n propagating 
in z direction, an electric field E™ polarized in x direction, 
and a magnetic field H ™ along the y direction. 
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The Mie theorjA 2 ^ of light scattering by a spherical 
core-shell particle has been worked out by Aden and 
Kerkei® long time ago. It can be found in many text- 
book^E. The geometry of the electromagnetic scatter¬ 
ing problem is shown in Fig. [l] An electromagnetic 
wave with wave number A~ 4 propagating in z direc¬ 
tion and an electric field polarized in x direction hits 
a particle centered in the origin of the coordinate sys¬ 
tem. The particle with total radius r 2 and refractive 
index iV 2 (A -1 ) = yje 2 (A~ 4 ) contains a core with radius 
7"i = fr 2 and refractive index iVi(A _1 ) = ^/£ 1 (A~ 1 ), where 
= E \ 2 (A _1 ) + * e i 2 (A -1 ) are fh e complex dielec¬ 
tric functions for the two regions and 0 < / < 1 is the 
filling factor. In the formulae below we use n = 27rA~ 4 
and w = 27 tcA - 1 instead of A -1 and v with c the speed of 
light, the abbreviations ki = nNi with i = l,2,m, where 
N m is the refractive index of the surrounding medium, 
and the size parameters a: 2 = 27rA~ 1 r 2 and X\ = 27rA _1 ri. 

In order to obtain the classical optical response of the 
particle we expand, as in the standard procedure, the 
electromagnetic fields outside and inside the particle in 
the vector spherical harmonics M, N and determine the 
expansion coeffic ients from the boundary conditions at 
the two interfaced 3 4 1? . Special attention is however paid 
to the shell region, where we follow XrP 

and Pena and 

Pal 24 and split the fields in a manner suitable to bring 
out the hybridization scenario. 


where Eq is the strength of the incident electric field. The 
radial dependence of these fields is described by Bessel 
functions of the first kind, indicated by the superscript 
(1), thus they can be interpreted as the penetrating fields 
(explaining the subscript p) in the surrounding medium 
(* = to). The scattered fields outside the particle may be 
written as 


ET = E E n (ia™N£l - b™M£l) , (4) 

n=1 

&T = — £ E n (ib™N 0 <£ + a™M£l) . (5) 

n=l 

The superscript (3) indicates, that Hankel functions of 
the first kind h n = j n + Wn are used for the radial de¬ 
pendence of the fields, whereby y n are Bessel functions 
of the second kind, which will be indicated by the su¬ 
perscript (2). Inside the particle the core and the shell 
region have to be distinguished. The penetrating fields 
inside the core (i = 1) are given by 




y(i)' 


A^c 


W/Ti n=1 
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and the fields inside the shell may be expressed as 

oo 

-Esheii = E E n (f n M° ln - ig n N^ n 

n =1 

+ V nMoll ~ iw nN^l) , ( 8 ) 

Hs hell = ~^t E n(9nM^ + if n N™ 
n = 1 

+ WnM^l + iv n N%l) . (9) 

The expansions ([8]) and (|9]) going back to Aden and 
Kerkei^ are the basis of most applications of the Mie 
theory to spherical co re-she ll particles. In our previous 
work we also used them 3 ®^. Their physical content how¬ 
ever is not obvious. To bring the physics to the forefront 
it is more appropriate to follow the hybridization sce- 
narijUHi ni a nd to express the fields inside the shell (i = 2) 
as follow^- 3 221 


Eshen = + Ep and H she ii = H 2 + Hp , (10) 

where A 2 , 7/ 2 and A 2 . are, respectively, the scattered 
fields (explaining the subscript s ) of a cavity consisting 
of core material and embedded in an homogeneous do¬ 
main of shell material and the penetrating fields of an 
homogeneous particle made out of shell material (hence 
the subscript p) embedded in the medium. Adopting 
Eqs. and Eqs. ([g]) ([ 7 ]) to the shell region, 


fourth multipole order (as are all coefficients in front of 
the spherical vector harmonics N e i). 

The hybridization scenario becomes apparent by writ¬ 
ing the expansion coefficients d 2 and a 2 in a form resem¬ 
bling the diagonal elements of a matrix resolvent 

G{E) -{ eA V- E ss-e) 1 ( 19 > 

describing the hybridization of two energy levels ea and 
Eb■ The basis for this matrix are the states |A) and \B). 
Lets consider 


Ga{E ) = (A\G(E)\A) 
1 


9a\E) 


\V\ 


( 20 ) 


9 b 1 (E) 


as a guide, where |E| 2 is the hybridization strength of 
the two levels and gi(E) = (E - £j) _1 is the resolvent of 
the isolated level i = A, B. The pole of gi{E) gives the 
energy of the noninteracting level i, while the two poles 
of Ga(E ) are the excitation energies of the interacting 
system. Once the coefficients are in this form their in¬ 
terpretation in terms of an hybridization scenario is thus 
obvious. 

Indeed the coefficients d 2 and a 2 obtained from the 


boundary conditions (15) (18) can be written as 



El=Y,En(ialN£l-blM£i) , 

n= 1 

(11) 


4 2 = — E MiblN™ + a 2 M e (3 >) 

w M2 n =l 

(12) 

and 

K =tM c lM ( 0 {l-id 2 n N£ n ) , 

n= 1 

( 13 ) 


Hi = E + iclN^l) . 

UP2 n =l 

(14) 

Since the total number of expansion coeffi¬ 
cients is unchanged, remaining in total eight, 
an,b™,al,bl,cl,dl,cl, and d 3 , the shell-medium 
boundary conditions at r = r 2 , 


rx(E™ + E™-El-E 2 p ) = 0, 

(15) 


rx(H™+H™-Hl-Hl) = 0, 

(16) 

and the 

core-shell boundary conditions at r = r\, 



rx(E 2 s+ El-El) = 0, 

(17) 


rx(Hl + Hl-H 1 p) = 0 , 

(18) 


are sufficient for their determination. For the derivation 
of the hybridization scenario we focus on the coefficients 
d 2 and a 2 which are resonant for small particles up to 


d? = 


X l( x 2 ) 


E l{x 2 ) - 


A^Qn )Zl(x 2 ) 

E\(x 1 ) 


( 21 ) 


and 


a„ = 


X^( Xl )X^x 2 )/E^x 2 ) 


E l (a?i) - 


X^\ Xl )Zl(x 2 ) 

El{x 2 ) 


( 22 ) 


with El, El, and X£ pl Z 2 playing the role of g/, g B x , 
and |E| 2 , respectively. The function X 2 has no direct 
analogue. For better readability we postpone the defini¬ 
tions of these functions to the point where they enter the 
discussion of the physical content of Eqs. (21) and (22). 


A comparison of (21) and (22) with (20) suggests, now 


at the level of the Mie coefficients, that the optical re¬ 
sponse of a core-shell particle is the response of two hy¬ 
bridized subsystems, with individual resonances deter¬ 
mined by 


El(x 2 ) = N 2 ^ n {N m x 2 )ip n (N 2 x 2 ) 

- X m £ >n {N m x 2 )4>' n {N 2 x 2 ) = 0 , (23) 


El(x 1 ) = Ni^N^x^ntNxXi) 

- N 2 ^ n (N 2Xl )^ n (N lXl ) = 0 , (24) 
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respectively, and an hybridization strength given by 
X^ >1 (xi)Z^(x 2 ). We have chosen all components to be 
nonmagnetic (yLti )2>m = 1). Notice in Eqs. (23) and (24) 
the Ricatti-Bessel functions V’n(p) = PJn(p) and £ n (p7 = 
ph n {p ) as well as the implicit definitions of the functions 
E 2 and E 4 appearing in (21) and (22). The physical 


content of (23) and (24) can be deduced by comparing 


these equations with the denominator of the scattering 
coefficient describing an homogeneous particle 4 


a„ = 


N^’^NmX^n^Nix) - N m il) n (N m x)ip' n (Nix) 
Ni£,' n (N m x)il) n (Nix) - NMN m xW n (N lX ) 


(25) 

which has to vanish in the resonance case. Here, N\ 
and N m are the refractive indices of the particle and the 
surrounding medium, respectively, and x is the size pa¬ 
rameter of the particle. Clearly, Eq. (23) is the resonance 


condition of an homogeneous particle with size parame¬ 
ter X 2 and refractive index N 2 embedded in a medium 
specified by N m , while Eq. (24) is the resonance condi¬ 


tion of an homogeneous particle with size parameter X\ 
and refractive index N\ embedded in a medium charac¬ 
terized by N 2 . If the real part of the dielectric function of 
the core e\ > 0 and the real part of the dielectric function 
of the shell e' 2 < 0 it is however the embedding medium 
which supports the surface modes. It is thus appropriate 
to consider (241 as the resonance condition of a cavity 
filled with material described by N\ and embedded in a 
medium characterized by N 2 . 

So far we extracted from the Mie formulae the two 
elementary building blocks, a shell-embedded core cav¬ 
ity and a medium-embedded homogeneous shell particle, 
whose surface modes hybridize to make up the optical 
response of the core-shell particle. We now turn to the 
hybridization strength X^ 1 (x\) Z 2 (x 2 ). It is interesting 
to analyze under what conditions it vanishes and what 
this implies for the expansion coefficients. The two sub¬ 
systems are noninteracting when at least one of the fol¬ 
lowing conditions holds: 


X n Pl ( x i) = IVi'0n(-N 2 Xi)'0( l (7Via;i) 

“ N2 , 'l>n(N 2 X 1 )4’ n (N 1 X i) = 0 , 
Zn(x 2 ) = N 2 £' n (N rn x 2 )£ n (N 2 x 2 ) 

- N rn ^ n (N rn x 2 )^ n (N 2 x 2 ) = 0 . 


(26) 

(27) 


Condition (26), which implicitly defines X ^ pl , is satis¬ 


fied when the refractive indices of the core and the shell 
are the same, N\ = N 2 . It leads to o? n - 0 and 

,2 = X n(x2) 

n ~ El{x 2 ) 

_ N 2 Cn(NmX 2 )4>n(N rn x 2 ) - N 2 t n (N m x 2 )i/j' n (N m x 2 ) 
N 2 C n (N m x 2 )ip n (N 2 x 2 ) - N m ^ n (N m x 2 )ilj' n (N 2 x 2 ) 

= 4 , (28) 


is now identical with the shell. Note, the numerator in 
the second equality defines the function X 2 appearing in 
Eqs. pll) and ((221). Thus, due to the absence of the core¬ 


shell interface, the scattering fields in the shell region 
vanish, leading to a 2 n = 0, and the penetrating fields in¬ 
side the shell become equivalent to the penetrating fields 
of the core, signalled by d\ - d 4 , with expansion coeffi¬ 
cients describing the fields inside an homogeneous parti¬ 
cle with refractive index N 2 embedded in a medium with 
refractive index N m . The core-shell particle is in this 
limit reduced to a medium-embedded homogeneous shell 
particle. 


The second condition (27), which implicitly defines Z, 


2 

_ 715 

reduces the core-shell particle to a cavity embedded in an 
homogeneous medium. In this case the refractive indices 
of the medium and the shell have to be identical, N m = 
N 2 , leading to = 1 and 

2 _ Ni'ijj' n (N 2 xi)ip n (Nixi) - N 2 ip n (N 2 xi)ip r n (NiXi) 
X l^' n {N 2 X 1 )llj n {N 1 X 1 ) - N^niNtXiWniN-LX-i) 

(29) 

Hence, due to the missing medium-shell interface, the 
penetrating fields of the shell become equal to the in¬ 
coming fields, which are considered to be the penetrat¬ 
ing fields in the medium, as can be seen by comparing 
Eqs. (13) and (14) for d' 2 rl = 1 with Eqs. Q and (2]). The 
expansion coefficient a\ is in this limit attached to the 
scattered fields of a cavity. Indeed, looking at the scat¬ 
tering coefficient of an homogeneous spherd^ a™ given 
by Eq. (25) and substituting N m -*■ N 2 and x -*■ x\ 
makes (25) identical to (29). Hence, for N m = N 2 the 


core-shell particle is reduced to a core cavity embedded 
in a shell medium. 

It should be noted that in the standard expansion the 
coefficients f n ,g n ,v n , and w n reduce in the respective 
limits also to the Mie coefficients of a medium-embedded 
homogeneous shell particle and a shell-embedded core¬ 
cavity 33 . What we have shown by working with a^, 6^, 
c^, and d ^ and rewriting them in a particular manner 
is that the two limiting cases are actually two build¬ 
ing blocks whose resonances are always virtually present. 
They are encoded in the functions E 2 (x 2 ) and E^(x i). 
The interaction X^ l {xi)Z 2 (x 2 ), controlled by the ge¬ 
ometry and the material parameters, defines their life¬ 
times and makes the optical response of the core-shell 
particle given by the poles of and a 2 n unique, in full 
accordance with the hybridization scenaricpSHHO] 

Before presenting numerical results for finite hybridiza¬ 
tion strength let us consider the case where the hy¬ 
bridization is turned off by a vanishing filling factor 


model 

A T o ( cm 

£o 

Soo 

core 

300 

3 

2 

shell 

600 

20 

2 


where d 4 is the coefficient arising in the expansions |6]) Table I. Parameters for the model core-shell particle, 

and Q for the penetrating fields inside the core which 
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Figure 2. (Color online) Real parts of the dielectric func¬ 
tions specified by Eq. ( |30| ) and the parameters in Table [I] as 
a function of the wave number A -1 . 

/ = ri/r 2 = Xi/x^. Depending on how the limit / -»■ 0 is 
taken, the core-shell particle reduces again to either one 
of its building blocks, a medium-embedded homogeneous 
shell particle or a shell-embedded core-cavity. If the fill¬ 
ing factor vanishes because the core becomes smaller and 
smaller, that is, because x\ -*■ 0, while X 2 is fixed, it is 
X^ x {x\) which vanishes. As discussed above, the core¬ 
shell particle is then reduced to an homogeneous particle 
made out of shell material and embedded in a medium 
characterized by a refractive index N m . If the filling fac¬ 
tor vanishes however because the particle becomes larger 
and larger, while the core size is fixed, that is, for X\ 
fixed and X 2 -* it is Z^(X 2 ) which vanishes and re¬ 
duces, according to the discussion given in the previous 
paragraph, the core-shell particle to a cavity made out of 
core material embedded in an homogeneous surrounding 
made out of shell material. For finite filling factors both 
X^’ l {x\) and Z%(x 2 ) are finite and mix the optical re¬ 
sponse of the building blocks. This will be discussed in 
the next section. 

The re-organization of the Mie theory discussed in this 
section for the particular case of a spherical core-shell 
particle can be applied to other composite particles as 
well provided the symmetry is high enough to yield an¬ 
alytical expressions for the expansion coefficients of the 
electromagnetic fi elds. In the Appendix we show this for 
a stratified sphere 18 21 " 23 -. 


III. ILLUSTRATION 

To illustrate the classical optical response of a spherical 
core-shell particle we consider a particle with a dielectric 
core and a dielectric shell. Previously we proposed to use 
this type of particle in a gas discharge as an electric probe 
with optical read-out 34 . The idea, formulated initially 
for an homogeneous dielectric particle 3 -^, is to determine 
the electric field at the particle’s position in the discharge 


Figure 3. (Color online) Solid lines show the wave num¬ 
bers of the bonding (subscript A) and antibonding (sub¬ 
script B) dipole (riA,B = 1), quadrupole (ha,b = 2), and 
hexapole (ua,b = 3) resonance as a function of the filling fac¬ 
tor / = ri/r 2 for a core-shell particle with radius r 2 = 0.6 pm. 
The dotted and dashed lines denote, respectively, the posi¬ 
tions of the resonances of the medium-embedded homoge¬ 
neous shell particle and the shell-embedded core cavity, while 
bullets indicate the positions of the maxima in the extinction 
spectrum. 


from the balance of the forces acting on it and the charge- 
dependent shift of one of its extinction resonances. The 
shifts can be maximized by localizing inside the shell the 
elementary charges the particle acquired from the plasma 
using materials with negative electron affinity as a core 
and materials with positive electron affinity as a shell. 

The extinction spectra we calculated with the goal of 
employing them as a charge diagnostics in a plasma are 
based on the complex dielectric functions of the real ma¬ 
terials®^. In our work on topological aspects of light 
scattering by dielectric core-shell particles 33 ^ we found 
however dissipation to blur higher order bonding and an¬ 
tibonding resonances. For the present purpose we employ 
therefore dissipationless model dielectric functions with 
parameters chosen such that bonding and antibonding 
resonances can be clearly identified up to the third mul¬ 
tipole order. 

As in our work on topological aspects of light scattering 
by dielectric core-shell particles®*, we set N„, = 1, that 
is, embed the particle in vacuum, and use 

e ' - £00 + A t 2 0 -^, e" = 0 , (30) 

A to a 

with A^o, the wave number of the transverse optical 
phonon, and £q and £«,, the dielectric constants at large 
and small wave numbers, given in Table [T] for the core 
and the shell, respectively. The real parts of the model 
dielectric functions are shown in Fig. [2] Of particular 
interest is the spectral range where surface modes at the 
core-shell interface are excited, that is, the range of wave 
numbers where the real part of the dielectric function 
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n =1 


n =2 

c 


n c =3 


Figure 4. (Color online) Energetic positions of the reso¬ 
nances of the medium-embedded homogeneous shell particle 
(■ n 3 = 1,2,3), the shell-embedded core cavity (n c = 1,2,3), 
and the core-shell particle ( ua,b = 1,2,3) for / = 0.3 and 
T 2 = 0.6 pm. From the decreasing deviations of the hybridized 
from the non-interacting resonances can be inferred that the 
hybridization strength X^ 1 (xi)Z%(x 2 ) decreases with in¬ 
creasing multipole order. 


of the shell is negative and the real part of the dielec¬ 
tric function of the core is positive. For the model di¬ 
electric functions plotted in Fig. [2] this is the case for 
600 cm -1 < A -1 < 1900 cm -1 . 

The resonances in this spectral range are depicted in 
Fig. 0 as a function of the filling factor / = ri/r -2 for a 
particle with radius r 2 = 0.6 /im. In accordance wit h the 
hybridization scenaricf- 830 and our previous worlf 1 ' 1 3 £ 
we label the resonances at lower wave numbers bonding 
(subscript A) and the resonances at higher wave numbers 
antibonding (subscript B). Justification for this labeling 
comes from the polarities of the induced surface charges 
at the two interfaces which are in-phase for the bond¬ 
ing and out-of-phase for the antibonding resonances 33 . 
The solid lines give the positions for the bonding and 
antibonding dipole (tia.b = 1), quadrupole (tia,b = 2), 
and hexapole ( n a,b = 3) resonance as obtained from the 
poles of Eq. while the dotted and dashed lines give, 
respectively, the solutions of Eqs. (231 and (24), that is, 
the positions of the resonances of the medium-embedded 
homogeneous shell particle and the shell-embedded core 
cavity. Also shown are the positions of the maxima in the 
extinction efficiency, Q t = ^ T,n=i(^ n + l)Re[a™ + 
where a™ and b™ are the expansion coefficients of the 
scattered fields A™ and H™ (see Eqs. Q and ([5])). The 
resonances in the extinction spectrum depicted by the 
bullets follow precisely the resonances in the expansion 
coefficients dj) and a 3 indicating that the denominator of 
a™, which is the resonant coefficient in the expansion of 
the scattered fields, can be written in the same form as 
the denominators of d^ and a^. 

The main feature of Fig. [3] is the increasing splitting 
between the bonding and antibonding resonances with 
increasing filling factor / = ri/r^. As expected from 


the hybridization scenario the hybridization strength 
X'^ ,l {xi)Z^ l {x 2 ) controlling the splitting increases with / 
because the decrease in shell thickness leads to a stronger 
coupling between the surface modes at the core-shell and 
the shell-medium interface. For / r; 0.2 the resonances 
of the core-shell particle merge with the resonances of its 
building blocks, implying that the hybridization strength 
becomes vanishingly small. Notice, for tia.b < 2 the 
bonding resonances merge with the resonances of the 
medium-embedded homogeneous shell particle, that is, 
the solutions of E^{x 2 ) = 0, while the antibonding reso¬ 
nances merge with the resonances of the shell-embedded 
core cavity, that is, the solutions of E*(x i) = 0. For 
n-A,B - 3, however, it is the other way around. The rea¬ 
son is simply the energetic ordering of the resonances of 
the building blocks. For n\ ,b < 2 the resonances of the 
shell-embedded core cavity are energetically above the 
resonances of the medium-embedded homogeneous shell 
particle, while for tia.b = 3 the ordering is reversed. Fig¬ 
ure [4] displays this ordering for / = 0.3. 

Another way to visualize the hybridization scenario for 
the core-shell particle is to look at the electric fields inside 
the shell. This is shown in Fig. [5] for the dipole {n = 1), 
quadrupole (n = 2), and hexapole (n = 3) resonances of a 
particle with radius V 2 - 0.6 pm and filling factor / = 0.3. 
The lower and upper panels in Fig. [5] depict respectively 
the bonding and antibonding resonances. That the as¬ 
signment is correct can be inferred from the black arrows 
visualizing the orientation of the electric field inside and 
outside the particle. From them the polarity of the in¬ 
duced surface charges follows verifying that it is in-phase 
for the bonding and out-of-phase for the antibonding res¬ 
onances. The wave numbers are tuned to maximize the 
spatial extension of the scattered fields outside the parti¬ 
cle, E ™ and ta king as a measure the distance of the 
singular points^- 3 ' in the outer Poynting field from the 
center of the particle (for a discussion of singular points 
in the dipole fields of dielectric core-shell particles see 
RefP). From the ratio of the intensity of the scatter¬ 
ing electric field inside the shell \E %| 2 to the intensity of 
the overall electric field inside the shell \Ep + E ^| 2 , shown 
by the color coding, can be moreover deduced that the 
influence of the scattering fields is for the antibonding 
hexapole resonance smaller than for the bonding one, in 
contrast to the dipole and quadrupole resonances where 
it is reversed. The reason is again the energetic order¬ 
ing of the resonances plotted in Fig. |4j The hexapole 
resonance of the shell-embedded core cavity drops below 
the hexapole resonance of the medium-embedded homo¬ 
geneous shell particle. As a result, for n = 3 it is the 
bonding resonance which acquires more cavity character, 
and is thus more affected by the scattering fields, and not 
the antibonding one as it is the case for n= 1,2. 

For core-shell particles with dissipation, having 
complex dielectric functions, the higher order bonding 
and antibonding resonances may be blurred. For a 
Ca 0 /Al 2 03 particle, for instance, only the dipole 
and quadrupole resonances can be clearly identified® 3 . 
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Figure 5. (Color online) Electric field distributions in the *2 plane of a core-shell particle (r 2 = 0.6 //m, / = 0.3) near the 
bonding (lower panels) and antibonding (upper panels) dipole (n = 1), quadrupole (n = 2), and hexapole (n = 3) resonances. 
The wavenumbers are A -1 = 1313.62 cm -1 and A -1 = 1607.85 cm -1 for the bonding and antibonding dipole, A -1 = 1455.5 cm -1 
and A -1 = 1530.22cm -1 for the bonding and antibonding quadrupole, and A -1 = 1488.706cm -1 and A -1 = 1511.805cm -1 for the 
bonding and antibonding hexapole resonance. The direction of the total electric field is shown by black arrows and the color 
coding gives the intensity of the scattering fields in units of the intensity of the overall electric fields inside the shell. Black 
circles indicate the particle and its core. 


Hence, higher order bonding and antibonding resonances 
may be distinguishable only for a judicious choice of 
materials. The multipole where the character of the 
bonding and antibonding resonance changes depends 
also on the materials. It can occur already for the 
dipole 14 . The character of the resonances determines 
the number and locations of the singular points as 
well as the spatial distribution of the dissipation in¬ 
side the core-shell particle 43 . Changing it by tuning 
the resonances of its building blocks may thus be in¬ 
teresting from a physics and an application point of view. 


IV. CONCLUSIONS 

The hybridization scenario for the classical optical re¬ 
sponse of composite objects states that it can be under¬ 
stood in terms of interacting surface modes arising at the 
interfaces of the object. 


For the particular case of a spherical core-shell particle 
we derived the scenario directly from the expansion of the 
electromagnetic fields inside and outside the particle in 
terms of vector spherical harmonics. By re-organizing the 
fields inside the shell region into the penetrating fields of 
a shell particle and the scattered fields of a core cavity we 
were able to derive formulae for the expansion coefficients 
resembling the diagonal elements of a matrix resolvent 
describing two hybridized energy levels. From the formu¬ 
lae the building blocks of the hybridization scenario, the 
shell-embedded core cavity and the medium-embedded 
homogeneous shell particle, as well as the coupling be¬ 
tween them could be straightforwardly identified. The 
physical content of the coefficients became thus immedi¬ 
ately clear. In an Appendix we also demonstrated that 
the same strategy can be applied to a stratified sphere 
containing an arbitrary number of shells. Re-organizing 
the expansion coefficients for the fields yields again ex¬ 
pressions which resemble the diagonal elements of a ma¬ 
trix resolvent describing a chain of hybridized energy lev- 


























els. The basic building blocks of the hybridization sce¬ 
nario and their coupling could thus be identified directly 
from the expansion coefficients. 

Initially the hybridization scenario has been deduced in 
the electrostatic approximation for small spherical strati¬ 
fied particles without retardation using equations for the 
surface charges induced at the interfaces. The derivation 
we presented includes retardation and is thus applicable 
to spherical stratified particles of any size. In addition 
it provides a road-map for analytically deriving the hy¬ 
bridization scenario for other composite objects with high 
symmetry for which analytical solutions of the Mie the¬ 
ory are available. The physical mechanisms buried in 
these solutions may thus become apparent. 
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APPENDIX 



Hi 




E n 

n 


Figure 6. Geometry of light scattering by a stratified sphere 
with total radius r*, containing k layers with radii r, for 1 < 
i < k, whereby the core is counted as the first layer (i = 1). 
Refractive indices A; characterize each layer and the medium 
(i = k + 1 = to). The incident electromagnetic plane wave is 
described by a Poynting vector S™ propagating in z direction, 
an electric field E ™ polarized in x direction, and a magnetic 
field H™ along the y direction. 


In this Appendix we re-o rganiz e the Mie theory of 
a stratified spherical particl^SHIHlII t 0 obtain the hy¬ 
bridization scenario also for this situation. 

As indicated in Fig. [6] the particle contains k layers. 
The core is counted as the first layer i = 1, while the out¬ 
ermost layer and the surrounding medium are labelled 
i = k and i = k+1 = to, respectively. The label i is used for 
both the itli shell and its interface to the (i + l)th shell. 
The abbreviations of the main text are adopted. Each 
layer i is characterized by its radius r,;, its magnetic per¬ 


meability Hi and its refractive index 7Vj(A _1 ) = ^/^(A -1 ), 
where £i(A -1 ) = e'(A _1 ) +ie"(A _1 ) is the complex dielec¬ 
tric function for the respective region. The total radius 
of the particle is the radius of the outermost layer ry.. 
As before we use ki = kN z with i - l,...,k,k + 1, where 
Nk +1 = N m is the refractive index of the surrounding 
medium, and the size parameters xi = 2nX~ 1 ri. 

The expansions for the incident and the scattered fields 
outside the particle as well as for the penetrating fields 
inside the core, given respectively by Eqs. ([!]) and |2), 
Q and ([5]), and ([6]) and (IT]), remain the same, while 
the splitting of the fields (and their interpretation) in 
the shell of the core-shell particle, viz. Eq. (101, is now 
adopted to each shell 2 < i < k of the stratified sphere. 
Hence, 


K he n = K + Ei and = H* + W p , (31) 

whereby the scattered fields of the (i - l)th subsystem 


interact with the penetrating fields of the ith subsystem 


K* = V e (J (lit 1 ) _ mi 1 ) \ 

^p Zj rjn V n lu oln ta n Jv elnj > 


n= 1 


kjC 

UHi n=l 




f(l) 


+ ic l j 
^ tc rt Jv ol n) 


(34) 

(35) 


as schematically shown Fig. [7] A resonance arising at the 
*th interface, which-in the counting scheme introduced 
above-separates the ith from the (i + l)th shell, is thus 
described by a scattered field propagating into the (* + 
l)th shell and a penetrating field propagating into the ith 
shell. The overall field in the ith shell is a superposition 
of the scattering fields of the (i - l)th subsystem and the 
penetrating fields of the ith subsystem. 

Again the expansion coefficients are determined by 
the electromagnetic boundary conditions. For the shell- 
medium interface at r = r*,, 

r x (Ep +1 + 4 fe+1 - E k s - E p ) = 0 , (36) 

f x (H k+1 + H k+t - H k - H k ) = 0 , (37) 


while for the shell-shell interfaces at r = rt with 2 < i < k 


rx(Ei +1 + & p +1 = 0, (38) 

fx(Hi +1 + Hl +1 -Hi-W p ) = 0. (39) 




r(3) ’ 


he 


K = E 

n=l 


(32) 

(33) 


The boundary conditions for the core-shell interface at 
r = 7*1 remain Eqs. |l7| ) and (181. 

Before we identify the hybridization scenario from the 
expansion coefficients emerging from these boundary con¬ 
ditions we recall that the Mie approach works with the 
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Figure 7. Superposition of the electric fields in the shell 
region i. The resonance arising at the interface i between the 
ith and (i+ l)th shell is characterized by a scattered field El + 
propagating into the (*+ l)th shell and a field Ep penetrating 
into the ith shell. 


expansion coefficients whereas the hybridization scenario 
is an energetic picture focusing on energies (frequencies) 
for which the denominators of the expansion coefficients 
vanish. In deriving the energetic picture of the hybridiza¬ 
tion scenario from the expansion coefficients the question 
arises therefore on which coefficient one should focus. It 
turns out that any coefficient becoming resonant in the 


considered parameter regime can be taken. A fc-layered 
sphere containing k subsystems is described by the ex¬ 
pansion coefficients c l n , b 1 * 1 , d l n , and a^ +1 with i = 1 ,k. 
To extract the hybridization scenario we focus on d l n . It 
becomes resonant for small objects and is thus appropri¬ 
ate for our purpose. 

Our guiding principle is again the mathematical struc¬ 
ture of the diagonal elements of a matrix resolvent de¬ 
scribing hybridized energy levels. Instead of two energy 
levels, appropriate for the core-shell particle, we now have 
to consider a resolvent for k energy levels. Generalizing 
the approach described in Sect. [IT] we thus have to rewrite 
d l n in a form resembling the diagonal elements of the ma¬ 
trix resolvent 


G(E) 


f e 1 - E V 21 0 ... O' 

V-2i £2 ~ E I/32 : 

0 ••• 0 

: Vkk-1 

\ 0 ... 0 Vkk-1 £ k - E/ 


(40) 


describing a chain of energy levels with nearest neigh¬ 
bor interactions where i = 1, •••, k and Vjo = 14 +i,fc = 
0. Using renormalized perturbation theory they are given 
bj®! (|i) denotes the basis states) 


Gi(E) = (i\G(E)\z) = 


fh\E) 


\Vi 


i+li\ 


\Vi 


ii-1 1 


9ih(E) 


\Vi 


i+2i+l | 


\v k 


kk -1 


9I-\(E) 


\Vi-u-2 

|W 


21 


9l\E) 


9T\E) 


(41) 


where gi(E) = (E - £i) _1 is the resolvent of the ith non¬ 
interacting energy level. Its pole gives thus the energy 
of this level. The poles of Gi(E) on the other hand give 
the k excitation energies of the interacting system. Once 
the coefficient d l n is in this form its interpretation along 


the lines of the hybridization scenario is obvious as in the 
case of a core-shell particle. 

Choosing all components nonmagnetic, that is, setting 
Hi = 1 for i < k + 1, the coefficient d l n resulting from the 
boundary conditions can be indeed rewritten as 


<& = 


X^x k )Xt\x k . 1 y--Xl l (x l )/[E^x 1 ,---,x k )E^(x ir --,x k y--E, 




(ahr-Zfc)] 


E i(xi ) - 


X^(x 1 ,...,x. l )Z^\x i+1 ) 


X^ixu-.^x^Z^Xi) 


(42) 


, X^H Xl ,...,x i+ 1 )Z^(x i+2 ) , X^-\x l ,...,x i _ 2 )Z\~\x i . 1 ) 

E l n + 1 (x i+ 1)---- \ ry k / __ X K^Xi- 1)- 


X cpk ~\x l ,...,x k . 1 )Z k n {x k ) 

E %(x k ) 


X?Hxi )Z 2 n {x 2 ) 

K(* 1 ) 


that is, in a form closely resembling Eq. (41). Obviously 


Ey and X^ pl 1 Z l n play, respectively, the role of g i 1 and 
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|Vii_i| 2 . The functions Xl and in the numer¬ 

ator, having no direct analogue, are irrelevant for our 
purpose as they are not involved in the identification of 
resonance frequencies and coupling strengths. For better 
readability we again present the definitions of all these 
functions when they enter the discussion of the physical 
content of Eq. (42). 


The analogy between a chain of k hybridized energy 
levels and the optical response of a fc-layered sphere can 


be perhaps best seen by expressing Eq. (42) as a diagonal 
element of a tridiagonal matrix. Suppressing the size 
parameters Eq. (42) can be written as 


di = 


xlxl~'-x. 


E^Et 1,k -Ei +1 ’ i+2 ’- 


(AMI i) 


(43) 


with 


M = 


‘ El yjx^z* 0 _ 

yjx^zi El JxTzi 

0 ^X^Zl 


\[x^ kir Z k El ) 


(44) 


r 


Comparing (44) with (401 it becomes clear that the op- definition of A')( p! are 


tical response of a fc-layered particle can be interpreted 
in terms of k hybridized subsystems with individual res¬ 
onances determined by (z = 1 


El(xi) = Ni£ , n (N i+1 Xi)'ip n (N i Xi) 

- N i+ i£(N i+ iXi)ip' n (NiXi) = 0 


(45) 


and hybridization strengths given by 
^n Pl_1 ( a; i, • • • i x i-i)Z l n {xi). Notice in Eq. (|45|) the 


implicit definition of the function E l n (xi) which also 
appears in Eq. (|42|). Clearly, Eq. (45) is the resonance 


condition of an homogeneous particle with size param¬ 
eter Xi and refractive index TV, embedded in a medium 
specified by N i+ 

We now turn to the hybridization strength of the 
subsystem z, which in contrast to the core-shell par¬ 
ticle interacts with two adjacent subsystems. The 
interaction with the (i - l)th and the (z + l)th 
subsystem is given by A£ pl_1 (:ci,... , 2 q_ \)Z l n (xi) and 
X^ l {x\,... ,Xi)Zn 1 (xi+i), respectively. Here, 


KiXi) = Ni&iNMXiKniNiXi) 

- Ni+lZniNi+iXi^KNiXi) 


(46) 


for i = 2,..., k and 


X?{x 1 ,...,x i )=X*{x i ) 

. X2*- 1 (x 1 ,...,x^ 1 )U' n (x i ) 


(47) 


En (^-1 ; • • * j %i- 1 ) 

for i = 1,..., k - 1 with X )j p0 = 0. The functions in the 


Xn(Xi) = NilpKNi+iXi^niNiXi) 

~ Ei+i'ip n (Ni+-LXi')'lp n (NiXi') , 

U l n (Xi) = N i+1 li} n {N i+1 X i )£,' n {N i X i ) 

- Niip' n (Ni +1 Xi)£ n (NiXi) , 


(48) 

(49) 


and 


En ’ ( X 1 , . . . , Xi) — 

X^ i - 1 (x 1 ,...,x i - 1 )Zl(x i ) 


E l n (xi) + 


En "’ 1 1 (xi,...,x i -i) 


(50) 


where E l n {xi) is defined in Eq. (45). The function X £ pl 
depends on the size parameters aq,•••,Its superscript 
cpi, standing for composed, propagating and zth shell, 
distinguishes it from the function A pl which depends 
only on the size parameter aq. The recursion (50) gener¬ 


ates the second continued fraction in the denominator of 


Eq. (42). 


Of interest is also under what conditions the zth sub¬ 
system of the stratified sphere can be isolated and how 
this affects the expansion coefficient d l n . The decoupling 
is a two step process. In a first step, the interaction with 
the outer subsystems is turned off by choosing the outer 
shells and the embedding medium to be identical to the 
(z + l)th shell, that is, setting Nj = W+i for z+1 < j < k+ 1, 
which leads to Zl(xj) = 0 for i + 1 < j < k. In a second 
step, we turn off the interaction between the inner sub¬ 
systems by choosing the refractive indices of the inter¬ 
nal layers j for j < i and the refractive index of the zth 
layer to be the same. Hence, Nj = Ni for j < i, so that 
X^ l ^ 1 {xi ,... jZi-i) = 0. To analyze how the isolation of 
the zth subsystem affects the expansion coefficient d l n we 
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need also to look at the functions 


and 


X l n (xi) = Niip n (N i+1 Xi)£' n (N i+1 Xi ) 

- Nilp'niNi+xX^ZniNi+iXi) 


(51) 


El 


" k {x ll ... 1 x k ) = 


K(^) + 


E n \%1 j • ■ • 5 ) 


(52) 


entering the numerator in Eq. (|42j). The functions 
E l n (xi ), Zl(xj) and Xf pl (xi,..., a;,) are defined in 
Eqs. {45L |46| ), and ( |47| ), respectively. Note, the re¬ 
cursion (52) defines the first continued fraction in the 
denominator of J42~ 


Due to the decoupling, that is, the 
particular identification of refractive indices, the inter¬ 
action term X^‘ l {x\,.. ., Xi)Z^ 1 {xi+\) in Eq. (52) van¬ 
ishes leading to E%"’ k (xi,---,Xk) = E l n (xi). In addition 
X° n (xj)lE° n {xj) = 1 for j > i. Thus, Eq. (42) becomes 


di = 


Xfai) 

E i(Xi) 


(53) 


with Xl(xi) and E l n (xi) given by Eqs. (51) and (45). 


Due to the decoupling d l n should be attached to the 
penetrating fields inside an homogeneous particle charac¬ 
terized by Ni embedded in a medium with the refractive 
index fV i+1 . Indeed, looking at the penetrating coefficient 
of an homogeneous sphere given by Eq. (281 and substi¬ 
tuting N m -*■ N i+ 1 , N 2 -*■ Nt, and x 2 -»■ Xi makes (28) 
identical to ( |53| . Hence, if Nj = N t for all j < i and 
Nj - N i+1 for all j > i + 1 the fc-layered particle is reduced 
to an homogeneous particle with radius r,; embedded in 
a medium, whereby the used materials are described by 
Ni and 7V i+ i, respectively. 


As in the case of the core-shell particle, the original 
expansion coefficients of a fc-layered sphered reduce 
in the respective limits also to the coefficients of the sub¬ 
systems. Due to the rewriting we demonstrate however 
that the subsystems are always virtually present. The 
rewriting identifies thus the limiting subsystems as basic 
building blocks whose coupling yields the unique optical 
response of the composite object. 
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